1 


Comment on "Failure of the Work-Hamiltonian 
Connection for Pree-Energy Calculations" 

J. Horowitz and C. Jarzynski 

University of Maryland, College Park 20742 

Nonequilibrium work relations establish a connection 
between nonequilibrium work values and equilibrium free 
energies. In a recent Letter, Vilar and Rubi [l| (VR) 
argue that the definition of work used in these relations. 
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is incorrect, and therefore the relations themselves are 
fundamentally flawed. In our investigations, however, 
we have reached the opposite conclusion 0, as have Im- 
parato and Peliti in direct response to VR. 

In Eq. [Jl A represents generalized coordinates 
(fli, a2, • • • ) describing the external bodies that we ma- 
nipulate to act on the system of interest; W is the in- 
tegral of force, multiplied by the displacements of these 
bodies Referring to Refs. for a broader discus- 

sion, in this Comment we illustrate that W has exactly 
the properties we associate with thermodynamic work. 
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FIG. 1: The Hamiltonian functions and equilibrium distribu- 
tions (shaded) at / = and f = fo- 

We consider the model system analyzed in Ref. 
described by the Hamiltonian 


f) = -kx^ - fx, 


(2) 


and evolving under Langevin dynamics with friction co- 
efficient 7 = 1. The force f{t) is switched on uniformly. 


from /(O) = to /(r) = /o; outside the time interval < 
t < T the force is held constant. This system evolves from 
equilibrium state A in the distant past (/ = 0), to equi- 
librium state B in the distant future (/ — fo). The initial 
and flnal Hamiltonian functions and canonical distribu- 
tions are shown in Fig.[TJ Since H(x; f{t)) is constant for 
t < and t > T, its value during these times is identiflcd 
with the energy of the system [l| . Using the equilibrium 
distribution p oc exp(— /3i7), we compute the internal en- 
ergy (E = J pH) and the entropy (S = — J plnp) for 
states A and B: Ea = (2/3)-\ Eb = (2/3)-^ - /o/(2fc), 
Sa = Sb = [1 ~ ln(/3fc/27r)]/2. From the thermodynamic 
definition of free energy, G = E — ST , we then get 


AG = Gf 
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(3) 


The negative value of AG reflects a decrease in internal 
energy, with no change in entropy (see Fig. [1]). 

For this model, the distribution of work values over 
an ensemble of realizations of the process, p{W), can be 
obtained using the approach of Ref. [§] . This distribution 
is a Gaussian with mean and variance. 
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This result implies that: W AG for every realization 
in the reversible limit (r oo); (W) > AG in the ir- 
reversible case (finite t); and {e~^^) — e~^^'^ for any 
value of T. Thus the work defined by Eq. [7] is consistent 
with the second law of thermodynamics, and its exponen- 
tial average correctly gives AG, when the free energy is 
defined by the expression G ~ E ~ ST. 

By contrast, VR obtain AG > for this model (Eq. 4 
of Ref. (H), and assert that a negative value would be 
"inconsistent with a nonspontaneous process" . We dis- 
agree. An undisturbed system certainly seeks to minimize 
its free energy (e.g., after the removal of a constraint), 
but when an external agent varies a parameter of the sys- 
tem, such as the field / above, then there is no universal 
restriction on the sign of the free energy change. For in- 
stance, by manipulating a piston we can either increase 
or decrease the Helmholtz free energy of a gas, according 
to whether we compress or expand the gas. 
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